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Previous analysis of a Paris turbulence experiment [1,2] 
shows a transition at the Paylor Reynolds number Re\ ftj 700. 
Here correlation function data is analyzed which gives further 
evidence for this transition. It is seen in both the power spec- 
trum and in structure function measurements. Pwo possible 
explanations may be offered for this observed transition: that 
it is intrinsic to the turbulence flow in this closed box experi- 
ment or that it is an effect of a change in the flow around the 
anemometer. We particularly examine a pair of "probe ef- 
fects" . Phe first is a thermal boundary layer which does exist 
about the probe and does limit the probe response, partic- 
ularly at high frequencies. Arguments based on simulations 
of the response and upon observations of dissipation suggests 
that this effect is only crucial beyond Re\ ftj 2000. Phe sec- 
ond effect is produced by vortex shedding behind the probe. 
Phis has been seen to produce a large modification in some 
of the power spectra for large Re\. It might also complicate 
the interpretation of the experimental results. However, there 
seems to be a remaining range of data for Re\ < 1300 uncom- 
plicated by these effects, and which are thus suggestive of an 
intrinsic transition. 

PACS numbers: 47.27.-i, 47.80,+v 

(submitted to Physical Review E) 



I. INTRODUCTION 

Recent velocity measurements in highly turbulent he- 
lium gas flow by Tabeling et al. ("the Paris group") [1-5] 
reveal a transition in the turbulent behavior. Their re- 
sults show that beyond some crossover Taylor Reynolds 
number Re\ ps 700 the flatness of the velocity derivatives 
ceases to increase. Such a transitional behavior has never 
been reported in open flows [6,7], while earlier results on 
closed flows, also in helium [8-10], have also revealed the 
existence of a transition at large Reynolds numbers. For 
the Paris experiment, it has been argued [5,11] that the 
crossover signals the instability of vortex tubes (worms) 
[5]; one may speculate also whether it signals the onset 
of K41 turbulence [12]. 

This paper is devoted to a further analysis of this 
crossover, including two main questions: 1. Can the 
transition be seen in the further analysis of correlation 
functions? 2. Can it possibly be understood as some 
kind of effect of the probe? 

To discuss the former issue we analyze structure func- 
tions D n {r) = ((u(x + r) — u(x)) n ) of order n=2,4,6 and 
the energy spectrum. We try to superpose structure func- 
tions and spectra for different Reynolds numbers Re. In- 
deed, the superpositions show that there are different 
forms of the curves above and below the transition, but 
within each subregion the forms of the curves seem quite 
similar. The shifts necessary to superpose such curves 
give a measurement of characteristic length or frequency 
scales in each measurement. These scales, and their de- 
pendence upon Reynolds number, are important because 
they can help to give insight into the physical origin of 
the transition. 

The central point of this paper is to examine the pos- 
sibility that the apparent transition might be somehow a 
reflection of measurement imperfections in the probe. We 
do this by studying two different mechanisms for probe 
effects in some detail. 

There are two likely sources of difficulties in the mea- 
surement, one being that the probe might have insuffi- 
cient time resolution, the other being that it might be 
too large. To study the first possibility, we follow the 
analysis of Grossmann and Lohse [14], who showed the 
possibility of loss of probe sensitivity in temperature mea- 
surements for Rayleigh Benard flow [15,8]. This probe 
effect occurs as a result of delays and averaging related 
to thermal diffusion through the partially stagnant gas 
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about the probe. In the Paris experiment considered in 
the present paper, the probe is a hot wire anemometer, 
which works by heating the gas around it. Once again, 
the probe measures temperature changes and can possi- 
bly have poor temporal response. This is the first probe 
effect we analyze in this paper. The second kind of effect 
we consider involves what happens when the probe size 
and the dissipative scale becomes comparable. One effect 
which can arise is vortex shedding behind the probe. In 
"ordinary" hot wire anemometry (i.e. cylindrical probes 
working in air), the operating conditions are such that, 
even at large (overall) Reynolds numbers, there is no vor- 
tex shedding behind the sensors [16]. In the Paris exper- 
iment, the Reynolds number based upon probe size is 
large enough so that it is reasonable to ask about vor- 
tex shedding. This issue has been partially addressed 
in the past [2,4]; the general conclusion is that if the 
vortex shedding frequency is outside the turbulent spec- 
trum, no big perturbation is expected. In the opposite 
case, the vortex shedding mode may couple with the tur- 
bulent fluctuations, and perturb the measurement. The 
vortex shedding is only one of several things which can 
happen when the probe has a size comparable with the 
dissipative scale. 

In outline then, in the next section, we show evidence 
for a change in behavior at high Reynolds numbers. The 
following section is devoted to looking at probe effects 
which might partially explain the observed change in be- 
havior. In the section 4, additional evidence is drawn 
from the measured Reynolds number dependence of char- 
acteristic scales. In section 5 we demonstrate by a sim- 
ple estimate that ISR quantities as scaling corrections or 
PDFs of velocity differences are hardly effected by the 
transition. The last section is devoted to conclusions. 

II. EVIDENCE FOR TRANSITION 
A. The Flatness 

The Paris turbulent velocity measurements are done 
in low temperature helium gas, following the idea of the 
Chicago convection experiment [17,15,8] which were built 
followed upon the experiments of Threifall [18]. The flow 
is driven by two counterrotating disks of radius R. Two 
different cells are used: Cell 1 with R = 3.2cm [3,1,2] 
and cell 2 with R = 10cm [1,2]. The velocity anemome- 
ter ("probe") of size d = 7/um is placed far enough from 
the boundary layers. The Reynolds number is defined as 
Re = Sli? 2 jv. The angular velocity £1 of the disks, which 
is about l-10Hz, remains about the same for all mea- 
surements whereas Re is varied by changing the helium 
pressure and thus the viscosity v. 

The transition was first observed [2] in the properties 
of the flatness F and the skewness of the velocity deriva- 
tive. (We focus upon the flatness since the transitional 



phenomenon was less visible in the skewness.). Fig 1 
shows a series of flatness measurements for the small and 
large cells, plotted as a function of Re\. This Reynolds 
number is experimentally defined in terms of the RMS 
fluctuations in the velocity. [1]. Fig 1 incorporates un- 
published data recently obtained with the smaller cell. A 
rough estimate for the experimental uncertainty in flat- 
ness is ±15%. At lower Reynolds number, the flatness 
F increases with increasing Re\. Then at a transition 
value of the Reynolds number, it seems to reach a peak 
and then level off at higher Re\. The peak flatness is 
located at Re\ comprised between 550 and 750 for the 
small and large cells (see fig. 1 and also fig. 5 of [2]). 
These peak values are identified as the Re\ number for 
some kind of crossover or transition. 

The nature of the transition can be seen by studying its 
characteristic length scale. In particular, the small scale 
length called li„f (defined in reference [1]) describes the 
boundary between the ISR and VSR as observed in the 
third order structure factor. The classical expectation 
is that a length like this will decrease as Re~ 3 l A . This 
length shows the expected behavior below a characteris- 
tic Reynolds number and instead it seems to saturation 
at higher Reynolds numbers. (See figure 5 in reference 
[1].) It was also found (see [1]) that the scaling range 
of apparent ISR behavior seems to saturate and maybe 
even shrink beyond the transition Reynolds number. 

B. Dissipation spectra and structure functions 

The previous work mostly focused on the scalings of 
the viscous dissipation. Here we seek further evidence 
of the transition and its nature. We first analyze the 
energy dissipation spectra I(k) = k 2 E(k) as a function 
of Reynolds number. (Here, E(k) is the usual power 
spectrum, obtained as a function of frequency and then 
translated into a dependence upon wave vectors through 
its mean velocity following the Taylor hypothesis.) We 
find that the curves for the dissipation spectra fall into 
two groups: There is one shape which holds for all curves 
below the transition and another shape which governs the 
curves above the transition. 

The spectra are calculated over 16 million points, cor- 
responding to 100-1000 large eddy turnover times for 
each graph. The data is then divided into about 2 9 se- 
quencies of 2 15 points. The power spectra are computed 
for each sequence and averaged over different segments. 
A Hanning window is used in computing the spectra [19]. 
To collapse the spectra, we determine the peak of the 
dissipation I p and its corresponding frequency k p follow- 
ing the approach of reference [20]. We then rescale the 
wavenumber and the dissipation spectrum by k p , and I p , 
respectively. The rescaled spectra are then plotted in 
Figure 2. The Figure 2a shows the collapse for different 
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Reynolds number below the transition, whereas the Fig- 
ure 2b for those above the transition. Figure 2c shows 
a comparison of the shapes for the two scaling curves. 
The solid line corresponds to a general fitting form taken 
from [20,13]: 



I(k) = I(k p ) 



ftp 



5/3 



(1) 



where the fitting parameters are a = 0.7,/? = —1, and /j, 
is determined to make the peak occur at k/k p = 1. In 
this fit, the second summand reflects the bottleneck en- 
ergy pileup at the borderline between inertial subrange 
and viscous subrange [21]. Note the good agreement be- 
tween the fitting formula and the data holds above the 
transition, but not below. Figure 2 then shows a sit- 
uation in which the dissipation spectra changes at the 
transitional Reynolds number, but has a rather constant 
form above and below. 

The relative shifts in the logarithmic scales define un- 
ambiguously the length and energy scale in the system. 
In a section 4, we shall study the Reynolds number de- 
pendence of the length scales in comparison with different 
theories of what might cause the apparent transition. 

To further test the existence of two groups of scal- 
ing functions, we examine different structure functions 
D p (r) = ((u(x + r) — u(x)) p ) in the same spirit. First 
we look at Z?2- Since D2 and the power spectra contain 
the same information, the results on D2 should provide 
us a double check on the previous finding. To collapse the 
different curves of D2 versus r, we use the same method 
as applied to the dissipation spectra. We find a charac- 
teristic distance for each plot by finding the point of the 
log-log plot at which the slope is a predetermined con- 
stant. The constant is picked so that the slope stands at 
a place in which the viscous dissipation is beginning to 
occur. We pick the slope equal to 1.5. (Note that the 
slopes in D2 vary from 2 to ~ 0.7.) We then shift the 
plots on the log scales by the coordinates of this fitting 
point and overlay the plots. The result of this analysis 
is shown in figure 3. Once again, we find the data sep- 
arate into two groups: one for Re\ < 800, and another 
for Re\ > 1500. The relative shifts in horizontal and 
vertical directions define the scale in length and velocity. 
We return to the length scale in section 4. 

We apply the same method to study D4 and Dq. We 
again collapse by fitting that at specific values of the 
slope in the log-log plot. In these cases, the chosen slopes 
are respectively 1.5 and 2.5. In figure 4, we see clear 
separation of two groups. And once again, the separation 
occurs between Re\ = 800 and Re\ = 1500. However 
the figure 5 for Dq do not show such separations. The 
failure of seeing two groups in Dq somewhat weakens the 
argument for a simple transition. 



Thus we have shown some additional evidence for the 
existence of a transition in the Paris experiment. But 
only D2, its Fourier Transform E(k), and D4 indicate a 
simple transition. The other correlation functions sug- 
gest a broadened transition. 

It is clear that there is some kind of change or transi- 
tion centered at Re\ of about 700. The question is, How 
does one explain the observed transition? Can it possibly 
be an effect of some behavior localized about the probe? 
We turn to those issues in the next section. 



III. PROBE EFFECTS 

In this section, we shall obtain a variety of order of 
magnitude estimates. To make these estimates, we will 
have to compare a characteristic dissipative frequency in 
the system to the characteristic inverse times produced 
by the probe itself. We estimate the dissipative frequency 
as u>d = U/(l0rj) in agreement with Zocchi et al.'s spec- 
tral measurements. This frequency is then the inverse of 
the time it takes for a disturbance of size 10?y to move 
past the probe. Here, we use classically expected [22] re- 
lation between the Kolmogorov length rj and the dissipa- 
tion rate e. The definition of the length is r\ = i^ 3 / 4 /e 1 / 4 
with rj having the approximate value 



1 = 30i?e" 3 / 4 , 
R 



(2) 



In our last steps of analysis we shall express our results 
in terms of the Taylor Reynolds number [1] 



Re x = 1.57Re 1/2 . 



From eq. (2) we obtain 



u d = — = ^ — Re 7 '\ 
10?y R 2 300 



(3) 



(4) 



Here and below, we do our estimates by setting the 
Prandtl number equal to unity. 



A. The thermal boundary layer 

Note that the results on li„f resembles those deduced 
from the temperature measurements in highly turbulent 
Rayleigh Benard flow [8]: Beyond a critical Rayleigh 
number Ra k, 10 11 , the scaling range of the tempera- 
ture power spectrum becomes smaller for increasing Ra 
and the measured dissipative power shows a weaker Ra 
increase as below the transition. Grossmann and Lohse 
[14] suggested that this apparent transition might really 
be an effect caused by the probe used to measure tem- 
peratures. The thickness 6 of the boundary layer around 
the probe sets a diffusive time scale 
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uj- 1 = 8 2 / K (5) 

where k is the heat diffusivity. Beyond u>s , the measured 
spectral strength will be reduced. If u>s is smaller than 
the UV spectral cutoff w^, this will affect the UV side of 
the inertial range in between u>s and wj. This explanation 
of the observed high Rayleigh number data as a possible 
probe effect is still unproven, but it has certainly never 
been ruled out. 

One might wonder whether a similar effect would af- 
fect the results of the Paris experiment. This specula- 
tion might be enhanced because simplest estimate for 
the value of Re\ at the transition of the Chicago exper- 
iment is (lO 11 ) 1 / 4 ?s 700. By some accident (?!) this is 
the same number as at the observed Paris crossover. 

Why is heat diffusion relevant for the velocity mea- 
surements? The anemometer is heated by an electrical 
current. The faster the fluid is passing by the probe, the 
more power is needed to keep the probe at some constant 
temperature, which is larger than the temperature of the 
surrounding. The probe gauge curve power vs. velocity 
is given by King's law and is experimentally known [3]. 

What velocities does the probe measure? It can not be 
the velocity directly at the probe as there will be a vis- 
cous boundary layer of thickness 8 around it in which the 
velocity is very small. The heat generated in the probe 
has to diffusively penetrate this layer. Thus the probe 
measures the velocity of the helium which is a distance 8 
away from the probe. The length scale 8 again sets the 
diffusive time scale described by equation (5). 

Our problem is to now estimate the important value(s) 
of 8 and then to see whether the u>s thereby generated 
provides an important cutoff on the responsiveness of the 
probe. 

This has been carried out in two different ways. One 
of us (VE) has done numerical simulations of the flow 
past the probe, assuming a laminar time-independent 
flow [23]. The calculation must be done numerically since 
the Reynolds numbers of the probe 

Ud d 

Reprobe — — R e ~^ (6) 

v K 

is of the order of 40 at the observed transition point. 
Here, d = l\im is the size of the probe. A steady state is 
reached in which the flow is constant and in which there 
is a constant flow of heat away from the probe. The cal- 
culated flow is then perturbed with a sudden upstream 
rise in temperature. The probe responds and the impor- 
tant result is that the probe response time is shorter than 
any time resolved in the experiment. The conclusion is 
that one can assume that the thermal response of the 
probe is perfect. 

The next question is: How large is the viscous bound- 
ary layer 81 There are, in fact, two answers. If the probe 
Reynolds number is of order ten to one hundred, then 
there are regions with a thin viscous boundary layer in 



which the flow comes very close to the cylinder. Ac- 
cording to Blasius theory this boundary layer thickness 
should be: 

8 oc d/y / Re probe . (7) 

The flow behind the probe produces a much larger region 
of stagnant fluid at rest, with a thickness comparable to 
the cell size. Thus we have also 

8 ocd. (8) 

The frequency produced by these lengths via equation (5) 
will only matter if they are smaller than or of the order 
of the Kolmogorov cutoff of equation (4). Since these fre- 
quencies differ by a factor of the probe Reynolds number 
they are quite different. The conditions they generate are 
different also. If the length d matters, then the frequency 
cutoff will bother us whenever the Reynolds number, Re 
obeys 

Re> (300 (£) J , (9) 

Conversely if the size of the Blasius boundary layer mat- 
ters, we will get the less stringent condition 

Re > (300i?/rf) 4/3 . (10) 

for the Reynolds number at which the thermal boundary 
layer of the probe becomes important. 

These answers are quite different. If the important 
distance is d, we get from equation (9) an estimate of the 
critical Taylor Reynolds number, Re\, as ps 1000 for the 
small cell and ps 2000 for the large cell. The experiment 
measures a crossover Re\ of order 700. These results are 
perilously close. 

On the other hand, if it is appropriate to use the Bla- 
sius length then equation (10) gives a much larger critical 
Reynolds numbers, in fact a critical value of Re\ which 
is greater than 10 4 and which can not be realized in the 
experiment. 

So which length should we use? To see the answer no- 
tice that the experiment measures the heat flow out of the 
probe. This heat flow is much larger in regions in which 
the boundary layer is thinner. Thus, for larger Re pro i, e 
the important regions are the ones which have the Blasius 
effect and are thinned by a factor of (iJeproje) -1 / 2 . The 
simulation [23] fully supports this point of view. Note 
also that this line of argument does not apply to the 
Chicago experiment [15,8] and its analysis [14]. In the 
Chicago case, the probe measures the average tempera- 
ture in its environment. Regions of thick boundary layer 
are as important as thin ones. Therefore, relatively slug- 
gish regions can effect the outcome. In contrast, in the 
Paris experiment, only the most responsive regions mat- 
ter. 
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B. The dissipation rate test 

One can construct a direct test of the probe respon- 
siveness. The energy dissipation rate - denoted by e - 
can be measured in two different ways [f]. One involves 
a measurement in the inertial subrange and determines 
a quantity which we call £isr- The other is a dissipative 
measurement and determines £visr- Both quantities can 
be expressed in dimensionless form by writing c as the ra- 
tio of e to U 3 1 R where U = Sli? is the large scale velocity. 
In this way one finds the two different (Reynolds number 
dependent [24]) ratios 

c € ,isr = tiSRjj^ and c € y S R = cvsrj^ (11) 

One can find cisr by using the fact that the third order 
structure function D^(r) obeys the Kolmogorov struc- 
ture equation [22] 

£>( 3 )(r) = -4er/5 (12) 

When the Paris group carries out their ISR measurement 
their results agree with the expected spatial scaling and 
thus enable them to construct c c jsr- This is plotted 
in figure 6. In addition, the energy dissipation rate e is 
measured by a method which uses the viscous subrange, 
£vsr, namely by determining e from the spectrum E(k) 
as 

/>oo 

e VSR = lf>v / dkk 2 E(k). (13) 
Jo 

Here, full isotropization has been assumed. Strictly 
speaking, cvsr is only based on ((9i«i) 2 ) which is (via 
Taylor's hypothesis [25]) the only experimentally accessi- 
ble contribution to the strain tensor diUj. This determi- 
nation of e then gives the other dimensionless constant, 
Cc,vsr, which is also plotted in figure 6. If the probe 
response was cut off at high frequencies, one would ex- 
pect c e ysR to be substantially smaller than c c jsr- The 
data shows no support for that hypothesis. So we must 
conclude that, within experimental error (which is rather 
substantial), the high frequency response of the probe is 
satisfactory. 

C. Vortex shedding behind the probe 

In this section, we discuss the effect of finite probe size 
on spatial resolution. We focus particularly upon the 
effects of vortex shedding from the probe. 

Zocchi et al. [1] have observed a series of peaks at the 
high end of the power spectra. The exact origin of these 
peaks is unknown. Possible sources are: vortex shed- 
ding, vibrations of probe and its support. In more re- 
cent experiments the fiber has been strained at a ten- 
sile strength ten times larger than before and the peaks 



have mostly disappeared. This indicates that some the 
observed peaks were a vibration of the fiber. Here, we 
study the vortex shedding to estimate where it will occur 
and the effect it is likely to have. 

The frequencies of vortex shedding can be estimated 
as: 

fv = St 1 -^ (14) 

where St is the Strouhal number, which is typically 0.2 
[27]; Uo is the velocity of an ambient fluid, and d is the 
diameter of the obstacle. We focus upon the shedding 
from the probe, which has a diameter d ~ 7/um. 

Two conditions must be satisfied for the vortex shed- 
ding to be important. First it must be present. Vortex 
shedding appears when the probe Reynolds number is 
above ~ 40 [27]. In figure 8 we plot the probe Reynolds 
number Re pro i, e as a function of Re\. The horizontal 
dashed line corresponds the known [27] onset value 40. 
Its crossing with the experimental data gives a reading of 
transitional value in Re\, which is between 600 and 800 
for small cell, and 600-1100 for large cell. These ranges 
are close to the observed transition in flatness. Next, 
the shedding frequency must be in an observable range 
of frequency, by using typical values for Uo ~ lm/s in 
equation (14), we estimate the vortex shedding frequency 
to be about IbkHz, which is in the observable range of 
frequency. 

We then ask when does vortex shedding have an im- 
portant effect upon the flatness? Since the flatness is 
measured as a scale on the order of typically brj, we ex- 
pect the vortex shedding begin to affect the flatness when 
the vortex shedding frequency (eq 14) is comparable with 
the Kolmogorov frequency ( eq. 4). Under Taylor's hy- 
pothesis, this is equivalent to comparing the shedding 
wavelengthes l v and Kolmogorov length rj. The shedding 
wavelength is given by: 

'« = TT" ~ L2rf < ( 15 ) 

27T/t, 

In figure 9 we plot rj as a function of Re\ . Once again, we 
find that the line 1.2c? crosses the data around the same 
range of transitional value in Re\ as seen in the flatness. 
Thus, the vortex shedding shows promise of explaining 
the position of the transition. 

Next we want to know how does the vortex shedding 
affect the measurements of the flatness. Let us decom- 
pose the velocity signal U (t) as 

U(t) = V(t)+Asm(u v t) (16) 

where V(t) denotes the velocity in the absence of shed- 
ding and the other term might reflect, for example, a 
shaking of the probe with angular frequency u> v . We 
know that the intrinsic signal, V(t) gives a very large 
flatness, of order ten. On the other hand a sinusoidal 



5 



signal, like the term in A, will give a much smaller flat- 
ness. Thus, we should expect that such an additional 
term would tend to reduce the flatness as seen in the 
experiments. A careful calculation bears out this point 
[28]. 

However the figures 8 and 9 show some puzzling fea- 
tures. Note the probe Reynolds number is defined as in 
equation (6). For a fixed geometry we expect Re\ to be 
a fixed constant times Re 1 ! 2 . Thus all the black dots 
(corresponding to the small cell) in figure 8 should fall 
on one straight line with slope 1/2 and all the white dots 
should fall on another, lower by an amount correspond- 
ing to the log of the ratio of the cell sizes, log 3. So why 
do the large cell data fall onto two lines? Figure 1 and 
8 collects data measured in two distinct cells : one is 3 
cm in radius, and the other one is 10 cm in radius. They 
cover four series of experiments performed in the period 
1994-1995; there is no simple relation between Re\ and 
Re, because, for the same cell size, we may have sub- 
stantially different integral scales (due to the fact that 
the rims have not the same size from one series to the 
other); moreover, we have worked with different velocity 
fluctuation rates. All this explains why, for a given cell 
size, we have not a simple relation between Re\ and Re. 
However, if we fit the means of the two lines we get a 
reasonable value of the separation between the curves. 
This separation realizes our theoretical expectation that 
the two intersection points should have Re\ values which 
differ by a factor of the square root of the size ratio or 
roughly 1.7. 

On the other hand, according to figure 1 there is no ob- 
vious difference in the transition point between the larger 
cell and the smaller. This lack of difference would be 
expected if the transition were an intrinsic effect; it is 
not expected if it results from the probe size via having 
Rsprobe with a constant value at the transition. Thus, 
figure 8 seems to explain the position of the transition 
as a probe effect, but cannot satisfy us on the question 
of the size dependence. We do note that the uncertain- 
ties in the measurements are large, thus a possible size 
dependence maybe overlooked. 

Now we return to the question on the size of the probe 
relative to ry. The shedding wavelength (eq. 15) is inde- 
pendent of the Reynolds number. On the other hand, the 
maximum dissipative frequency gives a scale 10?y which 
decreases with Reynolds number as Re~ 3 l A . If 10?y be- 
comes comparable or smaller than d, the measurements 
of the dissipative quantities, such as the length measure- 
ments discussed in the next section, will be influenced by 
the finite probe size. Figure 9 shows that in fact the probe 
size is considerably smaller than 10?y so that it is likely 
that the vortex shedding will give such high frequencies 
that it is unlikely to affect the measured length or even 
the flatness(The reader will recognize that we are getting 
onto dangerous ground, since we are trying to distinguish 
between ry and 10?y in an order of magnitude argument). 



What does one expect to see in the length measure- 
ments if probe effect is relevant? As we have already seen, 
the vortex sheddings introduce to the system an external 
length scale which is independent of Reynolds number. 
This length scale becomes observable if the measurements 
of interest is on the order of the probe size. This suggests 
a saturation of li„f above the transition, which is consis- 
tent with Figure (5) in reference [1]. Similarly, it predicts 
that the length scales divided by the Kolmogorov scale 
increase with Reynolds number as Re\ 3 ^ 2 . The compar- 
isons will be done in the next section. 

Finally we like to remark that a sharper conclusion 
concerning the vortex shedding (or other finite probe size 
effects) can be drawn if we compare two experiments with 
and without vortex shedding at the same Re\. In other 
words, we want to find out whether the probe Reynolds 
number becomes a relevant parameter above the transi- 
tion. One hopes that at Re\ greater than the observed 
transition, by either varying the cell size or the probe size, 
we can change the probe Reynolds number such that it 
is much greater than the critical value 40 for one case 
and much smaller than the critical value for another. If 
the transition is independent of the two situations, then 
we can believe Re pro i, e is not a parameter determining 
the transition, and thus the transition is very likely to 
be intrinsic. So far the current available data , as shown 
in figure 8, do not yet provide us with such an ideal sit- 
uation. We see that for Re\ > 600, Re pro i, e is roughly 
the same for both systems, thus we are unable to distin- 
guish the two cases. We hope that future experiments 
will eventually clarify these ambiguities. 

IV. CHARACTERISTIC LENGTHS: 
MEASUREMENTS VS. PREDICTIONS 

In this section, we discuss the the characteristic lengths 
displayed by the turbulence data. Below the transition, 
we expect that the Kolmogorov dissipation length pro- 
vide the characteristic scale for all short-distance phe- 
nomena. Beyond the transition, we expect that another 
length might enter the problem. In both the vortex shed- 
ding approach and the thermal boundary layer (TBL) 
approach, this other length turns out to be of the order 
of the size of the probe 1 . In this section, we measure all 
lengths in units of the dissipation length. Therefore, we 
expect to see a constant value for this ratio below the 
transition. If the probe dominates the behavior by either 
of the mechanisms discussed here, we expect the ratio of 



J This result is obtained for the TBL by substituting equation 
(7) into equation (5) and multiplying by the integral scale 
velocity. The result is a length which is essentially the probe 
size. 
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characteristic length to dissipation length to increase as 

Our first method defines the dissipation lengths by 
identifying the maximum energy dissipation scale k p as 
we did to collapse the dissipation spectra. Figure 10 
shows the inverse of k p normalized by k v as a function 
of Reynolds number. The overall behavior can be fit by 
a constant as shown by the horizontal line. If the probe 
affect the scales when Re\ > 700) we might expect the 
ratio to slope upward as shown. Either a constant or an 
upward slope is equally well supported by these data. 

Similarly, we can determine other lengths scales by 
identifying the relative shifts used to collapse D2, D4, 
and Dq as discussed in section 2. These lengthes nor- 
malized by r) vs. Re are shown in figure 11. Again, a 
horizontal line is draw across them, and compared with 
the prediction of the two probe scenarios. These data 
do not support the notion of an upward slope. Hence 
they suggest that the scales of order 30 J7 are not affected 
by probes unlike the flatness measurements, which cor- 
respond to scales of brj. 

Next, we define the cross-over length to be the length 
r m corresponding to the minima of the 2nd log-log deriva- 
tive of D'2- Here, r msn should capture the cross-over from 
a r 2 scaling behavior at small r to a r° 7 scaling in the 
inertial range. The Reynolds number dependence of r msn 
normalized by rj is shown in figure 12 in comparison with 
the probe prediction. No evidence for a probe effect is 
seen. 

The Reynolds number dependence of all these length 
ratios seems to be consistent with a constant or the 
probe-predictions within the scatter of data. 

Another related length scale li„f defined in refer- 
ence [1] did show a saturation beyond the transitional 
Reynolds number, this is a suggestive of the existence of 
an external scale in the system. The plot of li„f is the 
figure 5 in ref. [1] 



V. CAN THE TRANSITION BE SEEN IN THE 
ISR? 

The transition in the flatness of d\Ui at Re\ c = 650 is 
rather pronounced, as seen from figure 1. How does this 
transition towards K41 in a VSR quantity reflect in pure 
ISR quantities? In this section we will see that for ex- 
perimentally reachable Reynolds numbers the transition 
can hardly be expected to be visible in ISR quantities as 
velocity structure function exponents and PDFs of ISR 
quantities. This finding is in agreement with the mea- 
surements of Tabeling et al. [4,2,1]. Thus the ISR results 
of the Paris experiment do not contradict the existance 
of the transition towards K41. 

For our estimate here it is sufficient to sketch the flat- 
ness F(r) = (v?)/(v 2 ) 2 as follows [29]: 



F(r) 



Fsat - 
— F 

J- OO 



: 10 

-2C2 



for 
for 
3 for 



r < I0r) 
10r] < r < L 
r>L 



(17) 



L is the integral length scale. From [1] we have L = 
4cm = 0.4i? for the large cell, independent of Re\. From 
the sketch (17) we immediately obtain the Re\ depen- 
dence of the scaling corrections to K41, 

8(Re x ) = \URex) - 2( 2 (Re x )\ = , ^^(o^ - ( 18 ) 



IgiL/W^Rex))- 



The ratio L/lOrj scales like 



L 

10?7 



cRe, 3 ' 2 . 



(19) 



Rather than taking c from eqs. (1) and (2) of section 
2 we adopt it to the experimental value of 8 for the 
Reynolds number Re\ = 650 of the transition. From 
[4] we have ( 4 = 1-25 and (2 = 0.70. Thus 6 = 0.15 and 
c = 0.185. The Re\ dependence of the scaling correc- 
tion 8(Re\) is very weak and approaches its K41 value 
only logarithmically, 8(Re\) oc l/lgi?eA, suggesting that 
l/logi?eA rather than l/Re\ is the small parameter in 
turbulence. Some numbers for still experimentally reach- 
able Reynolds numbers are given in table 1. This weak 
decrease is in agreement with the slight experimental de- 
crease of the scaling correction 8( n as shown in figure 8 
of [4]. 

Now we focus on the PDF of v r . In figure 2 of ref. 
[2] no detectable dependence of the PDF for v r for fixed 
scale r = 490//m on the Reynolds number was noticed. 
With our above sketch (17) of F(r) we readily calculate 



■ r \ -6 (Rex) 

3| l) 



(20) 



The Reynolds number dependence of F(r / L, Re\) for 
fixed r/L = 490 [im / 4cm = 0.012 again is very weak, 
table 1. 

Finally, we give as a characteristics of the PDF also 
the stretching exponent [3 in a parametrization 



PDF(v r ) oc exp 



(21) 



which is well known to fit the tails of experimental PDFs. 
[3 = 2 means Gaussian PDF. The flatness and the stretch- 
ing exponents are connected by [29] 



F 



r(i//?)r(5//?) 
(r(3//3)) 2 ' 



(22) 



The stretching exponent [3 again only very weakly de- 
pends on the Reynolds number, table 1. They must be di- 
rectly compared with figure 2 of ref. [2] where also hardly 
any dependence is seen. 
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To summarize this section: The Paris transition to- 
wards K41 can not be expected to be pronouncedly seen 
in pure ISR quantities as PDFs for fixed r or scaling 
exponents ( n . To look for a similar transition in other 
geometries one should thus focus on VSR quantities as 
e.g. the (hyper)flatnesses. 

VI. CONCLUSIONS 

The Paris experiments show a transition in behavior 
for Re\ pa 700. This transition can be seen in the flatness 
and indeed in measurements of individual moments like 
the power spectrum, D2, D4, and Dq. However, one 
might worry that the apparent transition was caused by 
the finite size or finite frequency resolution of the probe. 
We have analyzed herein two effects which may affect the 
measurement performed in the Paris experiment: a. the 
thermal response of the probe, which imposes a limitation 
on the temporal resolution and b. the vortex shedding 
behind the probe, which distorts the spatial resolution of 
the probe. What should we conclude? 

The two probe scenarios discussed in this paper both 
sound plausible. We cannot prove they are the direct 
cause of the transition, and in fact, we find partial evi- 
dence against each. In the end, we conclude that the mea- 
surements up to Re\ ~ 1300 are probably not perturbed 
by the two effects which are discussed in the present pa- 
per. 

We suspect the thermal boundary layer effect to be rel- 
evant at higher Reynolds number ~ 2000, but not in the 
region where the transitional behavior is observed. The 
key point here is the presence of a boundary layer which 
enhances the thermal exchanges between the fiber and 
the surrounding fluid. This is confirmed by numerical 
analysis [23], which showed that the thermal frequency 
response is outside the turbulent spectrum in a comfort- 
able range of Reynolds number around the transition. 
This is further checked by the observation that the dissi- 
pation rate seems correctly measured in the same range 
of Reynolds number; therefore the thermal response of 
the probe seems appropriate to the measurement of dis- 
sipative quantities. 

Vortex shedding probably do not affect the mea- 
surement in an appreciable range of Reynolds numbers 
around the transition, either. Although it gives the right 
transition value in Re\ and produces the observed de- 
crease of flatness with Re\, our basic argument against 
it is using the fact that the transitional Re\ for small and 
large cell are roughly the same within the experimental 
error. On the other hand, this argument is somewhat 
weakened due to the large error in the experimental mea- 
surements. We also remark that a set of improved exper- 
iments show less of the anomalous peaks in the spectra, 
and yet the transition persists. 



Perhaps neither the vortex shedding nor the thermal 
layer is in itself the right explanation. However, the tran- 
sition does occur when the probe Reynolds number is 
high (of order 40) and when the dissipation length is of 
the same order of magnitude as the probe size. There 
are likely to be many other possible effects, not explored 
here, which only depend on two essential ingredients: 1) 
the finite size of the probe, 2) the injection of energy in 
the small scale comparable to the size of the probe. Thus, 
another effect of finite probe size or response might in- 
tervene and produce a false signal of a transition. 

On the other hand it is entirely possible (and likely) 
that the transition observed is real and has nothing to 
do with probe effects. 

To more fully understand the nature of the transition, 
we will require further experiments with closed flows. 
They are likely to involve Helium and probes similar to 
those employed here. Clearly, it would be very desirable 
if a major piece of the next experiments were devoted to 
understanding further the behavior of the probes, and of 
the flow in their neighborhood, and how this flow evolved 
with Reynolds number. To fully realize the potential of 
the experimental method, we need further development 
of the technique for using probes like these. Perhaps one 
can probe the velocity and temperature field around large 
probes with tiny ones to further understand the probe ef- 
fects. 
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Table 1 

Reynolds number dependence of ISR characteristics for 
various Reynolds numbers beyond the transition. Only 



weak dependence is 


detected. 




Rex 


6 (Rex) 


F (490 fim, Rex) 


13 (Rex) 


650 


0.150 


5.82 


1.02 


1000 


0.139 


5.54 


1.06 


1500 


0.130 


5.32 


1.09 


3000 


0.117 


5.03 


1.14 


5000 


0.109 


4.85 


1.18 



FIG. 1. The flatness of the velocity derivatives is plotted 
against Rex- Results for the small cell are given as the black 
points and for the large one as the empty circles. This same 
convention is used in the later figures. Each point shown here 
is the average of three experiments performed at the same 
Rex, in the same cell, with the same probe. 

FIG. 2. Collapse of dissipation spectra. Parts a and b re- 
spectively show overlays of different dissipation spectra for 
Rex below and above the transitional value. Part c shows 
the comparison of two spectra, one from above and one from 
below. 

FIG. 3. Two groups of D2 curves. As in figure 2 this figure 
contrasts the relatively unvarying behavior above and below 
the transition with the somewhat larger change which occurs 
across the transition. In this and next two figures, we use 
open symbols to indicate the flows below transition, and solids 
above the transition, (chicago) 

FIG. 4. Two groups of D4 curves. As in figure 3 this figure 
contrasts the behavior within the region above and below the 
transition with the change which occurs across the transition, 
(chicago) 

FIG. 5. Trying to collapse of D% curves. As in figure 3 
this figure contrasts the behavior within the region above and 
below the transition with the change which occurs across the 
transition. This contrast is now constructed for correlation in 
the velocity cubed. However, we do not see the separation of 
groups as we did in the previous two figures, (chicago) 

FIG. 6. Dimensionless constants which measure the en- 
ergy dissipation rate, eps, found by two different experimen- 
tal methods. One method uses Integral subrange data, the 
other using viscous subrange data. The result is plotted as 
a function of the Taylor Reynolds number. Each point is 
the average of three experiments performed at the same rel, 
in the same cell, with the same probe. Notice that the two 
measurement agree within experimental error. Earlier mea- 
surements showed an apparent discrepancy in which the dissi- 
pative range quantity fell below the integral range one. Those 
measurements lead to the assumption that probe effects de- 
graded the data at high frequencies [26,14]. These plots show 
the Paris group's latest measurements performed at smaller 
fluctuation rates (around 20%, to be compared to 35% for [1]) 
, and do not reveal any discrepancy between eisR and evsR 

FIG. 7. The spectra for Rex = 1626. The arrows marks 
the positions of the estimated Kolmogorov wavenumber, and 
the of the vortex shedding. 
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FIG. 8. The various probe Reynolds numbers we achieve 
in the experiment for the two cells (small - black points - and 
large - circles -). The scatter is partly due to the fact that 
several fluctuation rates are considered in the compilation. 
The line shows the critical value, ~ 40, for the onset of vortex 
shedding. 

FIG. 9. The various ratios d/r] we achieve in the experi- 
ment. Same symbols and same remarks as in Figure 5a. 

FIG. 10. The relative shifts used to collapse dissipation 
spectra divided by k v vs. Re\. 

FIG. 11. The relative shifts used to collapse D2, Di and D% 
divided by r] vs. Re\. 



FIG. 12. 

' mm 

vs. Re\. 



[1] G. Zocchi, P. Tabeling, J. Maurer, and H. Willaime, 

Phys. Rev. E 50, 3693 (1994). 
[2] P. Tabeling et al, Phys. Rev. E53, 1613 (1996). 
[3] J. Maurer, P. Tabeling, and G. Zocchi, Europhys. Lett. 

26, 31 (1994). 

[4] F. Belin, P. Tabeling, and H. Willaime, Physica D (1996), 
in press. 

[5] F. Belin, P. Tabeling, and H. Willaime, J. de Physique 

(1996), in press. 
[6] C. W. van Atta and R. A. Antonia, Phys. Fluids 23, 252 

(1980). 

[7] For a comprehensive recent review see K. R. Sreenivasan 

and R. A. Antonia, "The phenomenology of small scale 

turbulence", Ann. Rev. of Fluid Mech., 1996, submitted. 
[8] I. Procaccia et al., Phys. Rev. A 44, 8091 (1991). 
[9] X. Z. Wu, L. Kadanoff, A. Libchaber, and M. Sano, Phys. 

Rev. Lett. 64, 2140 (1990). 
[10] B. Chabaud, A. Naert, J. Peinke, F. Chilla, B. Castaing, 

B. Hebral, Phys. Rev. Lett. 73, 3227 (1994). 
[11] J. Jimenez, A. Wray, P. Saffman, and R. Rogallo, J. Fluid 

Mech. 255, 65 (1993). 
[12] A. N. Kolmogorov, CR. Acad. Sci. USSR. 30, 299 (1941). 
[13] S. Grossmann and D. Lohse, Phys. Fluids 6, 611 (1994); 

Phys. Rev. E 50, 2784 (1994). 
[14] S. Grossmann and D. Lohse, Phys. Lett. A 173, 58 

(1993). 

[15] B. Castaing et al, J. Fluid Mech. 204, 1 (1989). 

[16] J. O. Hintze, Turbulence, Mc Graw Hill, 1975. 

[17] F. Heslot, B. Castaing, and A. Libchaber, Phys. Rev. A 
36, 5870 (1987). 

[18] D. C. Threifall, PhD. Thesis, Univ. of Oxford, 1976 (un- 
published); J. Fluid. Mech., 67, 1 (1975). 

[19] William H. Press, et. al., Numerical Recipes (Cambridge 
University Press, Cambridge, 1986) page 425. 

[20] Z. S. She and E. Jackson, Phys. Fluids A 5, 1526 (1993). 

[21] G. Falkovich, Phys. Fluids 6, 1411 (1994); D. Lohse and 
A. Miiller-Groeling, Phys. Rev. Lett. 74, 1747 (1995). 

[22] A. S. Monin and A. M. Yaglom, Statistical Fluid Me- 
chanics (The MIT Press, Cambridge, Massachusetts, 
1975). 

[23] V. Emsellem, C. R. Acad. Sci. Paris, in press (1996). 

[24] D. Lohse, Phys. Rev. Lett. 73, 3223 (1994); S. Gross- 
mann, Phys. Rev. E 51, 6275 (1995); K. R. Sreenivasan 
and G. Stolovitzky, Phys. Rev. E 52, 3242 

[25] G. I. Taylor, Proc. R. Soc. London A 164, 476 (1938). 

[26] K. R. Sreenivasan (unpublished). 

[27] D.J. Tritton Physical Fluid Dynamics (Oxford Scientific 

Publishers, second edition) page 28. 
[28] J. Wang, unpublished 

[29] S. Grossmann and D. Lohse, Europhys. Lett. 27, 347 
(1993). 



10 



